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Section 2: GENERAL APPROACH 


The new features uhat have been added to PEGASUS 5 are discussed in this section. 


2.1 Automatic Minimum Hole Cutting 

An automatic hole-cutting capability is provided with PEGASUS 5.1 In many cases, a good hole 
structure can be developed with no intervention by the user. However, the hole generation process can 
be manipulated through user inputs if necessary. 

The hole generation process requires a definition of the solid boundaries of the configuration. This 
information is supplied in the standard input file. See the $BCINP namelist for a description of this input 
and Appendix B for an example. This file can either be created manually using an editor, or can be 
generated automatically (with possibly some modifications necessary) if OVERFLOW input is 
available. 

The automatic hole generation process is illustrated using a two-dimensional example (i.e., a 3-element 
airfoil, used throughout this manual to illustrate aspects of PEGASUS operation). Figure 1 illustrates the 
solid boundaries of the airfoil, as defined in the boundary condition file. 







Figure 1. 1'hree- Element Airfoil Solid Boundaries 

The entire point of the hole-generation process is to partition the region surrounding the airfoil into 
"inside" and "outside" regions. In PEGASUS 5. 1 . this is accomplished using Cartesian meshes, where it 
is desired to mark each Cartesian element as an "inside", "outside", or "fringe" element. In the simplest 
case, a single Cartesian mesh is generated which contains the solid boundaries of the configuration. The 
elements of the Cartesian mesh which intercept the solid surface elements of the airfoil are found: Fthese 
are designated as fringe elements. Some of the fringe elements in the slat-wing region of the 3-element 
airfoil are depicted in Figure 2. 



Figure 2. Fringe Elements 

It is assumed that the comer elements of the Cartesian mesh are "outside" elements. Any unidentified 
(i.e„ non-fringe) element that is adjacent to an outside element must itself be an outside element. The 
outside region is then identified by marching from the comer elements inward until no more elements 
that are adjacent to outside elements can be found. The outside region is thereby completely defined, and 
is depicted in Figure 3. Finally, any element remaining that is not either an outside or fringe element 
must be an "inside" element. Inside elements are depicted in Figure 4. 




Figure 3. Outside Elements 



Figure 4. Inside Elements 

The Cartesian mesh is now a completed "hole map". Given an arbitrary grid point, the Cartesian element 
within the hole map in which the point resides can very quickly be identified. Points that are found to be 
encompassed by "outside" or "inside" elements are marked as field points or hole points, respectively. 
Points that fall within fringe elements must undergo further processing. PEGASUS 5. 1 uses a 
"line-of-sight" algorithm to determine the status of an arbitrary grid point: this algorithm simply states 
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element. The algorithm is illustrated in Figure 5. Points that can "see" an outside element are field 
points; points that can "see" an inside element are hole points. 
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Figure 5. Line-of-Sight Algorithm 

In the example of Figure 5, Point A is outside. Point B is inside, and will be marked as a hole point. 
Points that cannot be identified are considered outside points. 

2.2 Automation of Outer Boundary Specifications 

The automation of the specification of outer boundaries is relatively straightforward. If all boundary 
conditions have been specified for a mesh, the minimum and maximum index surfaces that have not 
been specified as boundary conditions for the flow solver are designated as the outer boundaries. Outer 
boundaries are assumed to receive boundary information via interpolation from other meshes. 

To automate the outer boundary specifications within PEGASUS, the boundary conditions for each grid 
must be available as input to PEGASUS. The boundary conditions are entered into PEGASUS using the 
SBCINP namelist in the standard input file. By default, PEGASUS will determine outer boundaries for 
all meshes. 

The manual approaches for specifying both outer boundaries and interior blanked regions are still 
acceptable inputs to PEGASUS. These inputs can still be used separately or in combination with the new 
method. 

In PEGASUS version 4, the user specifies whether a hole boundary is to have a single or double fringe 
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i. be v.vpiieiiiv specified. Now. with outer boundaries automatically specified, the fringe level for outer 
boundaries can be set with a single input. The user can set the fringe level for holes and outer boundaries 
with a single input, or set the fringe level for holes and outer boundaries separately. It should be noted 
that not all flow solvers can handle a mixed single/double fringe level of interpolated boundary points. 


2.3 Interpolation 

In version 4. the LINK list input is a priority list containing the names of meshes that are to be searched 
for interpolation stencils that can provide information for interpolated boundary points. The meshes in 
the priority list are searched in the order in which they appear in the input. If a legal interpolation stencil 
is found, then the processing for that interpolated boundary point is discontinued. If no legal 
interpolation stencil can be found in tnc first mesh in the priority list, then the second mesh is searched 
for an interpolation stencil. The process continues until either (1) interpolation stencils have been found 
for all interpolated boundary points, or (2) the mesh priority list is exhausted. Any interpolated boundary 
points that have no legal interpolation stencil in any of the meshes in the priority list are termed 
'orphans/' 


The user typically places the finer meshes at the top of the priority list and coarser meshes at the bottom 
of the list. This method usually gives adequate results but still leaves the user with only visualization 
techniques to evaluate the size difference between the interpolated boundary point cell and the 
interpolation stencil. (An interpolated boundary point cell is the cell adjacent to the interpolated 
boundary point in the positive J, K, and L directions.) It is assumed here that to obtain the "best" 
interpolation, the interpolated boundary point cell and interpolation stencil cell must have similar sizes. 
Use of the link-list approach, although simple to implement and computationally efficient, may not find 
the "best" pairing of an interpolated boundary point cell and an interpolation stencil. For example, if a 
user specifies a coarse mesh at or near the top of the link list, a significant difference in size may be 
found between the interpolated boundary point cell and the interpolation stencil. 

Automation of this process must therefore take into account the sizes of the interpolated boundary point 
cell and the interpolation stencil when determining the "best" pairing. Additionally, the user needs to 
receive diagnostic information to determine the worth of the interpolated boundary points-interpolation 
stencil pairs. 

The interpolation stencil that is determined to be the "best" for an interpolation boundary point is based 
on quality and the cell difference parameter. Quality is defined as a measure of the information that is 
being interpolated from non-hole points. As in previous versions of PEGASUS, if one of the stencil 
points is a hole point the interpolation stencil is not accepted. Quality is computed by first creating a 
temporary cell with the values at each of the eight points set equal to 1.0 (the point is a field point) or 0.0 
(the point is a boundary point). Quality is then determined by performing a tri-linear interpolation using 
the computed interpolation coefficients. If all of the stencil points are field points, quality will have a 
value of 1 .0. If for example, one of the stencil points is an interpolated boundary point and the 
interpolation coefficients are 0.5, 0.5, and 0.5, and, quality would be equal to 0.875. 

The cell difference parameter (CDP) is defined as follows: 
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Values for the cell difference parameter will vary from 0 (the best) to "big." It is yet to be determined 
what is "too big." Examples of the cell difference parameter are shown in Figure 6 for 2-dimensional 
cells. 



Figure 6. Examples of Cell-Difference Parameter 


Quality and the cell difference parameter are used to determine from which interpolation stencil, and 
therefore, from which mesh, an interpolated boundary point receives its data. Figure 7 depicts this logic. 
The quality value is checked first. If the quality is less than a user specified cutoff level (QCUTOFF), 
the interpolation stencil is no longer considered. 
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Figure 7. Logic for Interpolation Stencil Choice 


If the quality of the interpolation is better than a previously selected interpolation stencil by some 
tolerance (QTOL), the newer interpolation stencil is selected. If the value of quality is within the 
tolerance value, the interpolation stencil with the lesser value for cell difference is selected. If the 
computed quality is less than the current quality minus the tolerance, the previously accepted 
interpolation stencil is kept. 


2.4 Level 2 Interpolation 


Once the minimum hole has been established, the hole typically must be enlarged to improve the 
communication among the meshes. When only a minimum hole is cut, the fringe about the hole often 
causes very coarse mesh points to interpolate from very fine mesh elements. This disparity in mesh 
resolutions can severely impact accuracy, especially in viscous solutions, where a coarse mesh might be 
attempting to interpolate flow field information from within a boundary layer. In addition to being able 
to establish blanked regions in meshes with solid boundaries, it is also a requirement to have the 
capability to create holes with refinement meshes that are added to the domain. This cannot be 
accomplished through the usual hole-cutting methods, since refinement meshes typically do not have 
solid walls that would be used for automatic hole cutting. 


In order to handle these requirements a new approach for creating the interpolation boundary was 
developed. This approach, termed "level 2" interpolation, has the basic philosophy that coarser mesh 
points are interpolated from finer mesh points wherever possible. If more than two meshes are present, 
the coarser mesh will interpolate from the medium mesh, which will interpolate from a finer mesh. For 
all interpolations of this type quality is kept equal to 1 .0. This approach results in holes that are enlarged 
only if better communications are obtained by doing so. Therefore, the flow field will be computed for a 
given region only by the mesh with the best resolution. In addition, proper communications is 
established by this approach when refinement meshes are added. Figure S depicts the hole structure in 
the outer mesh in the 3-element airfoil case. Note that the outer boundary of the wing mesh has also 
been moved in so that the outer boundary of the wing mesn obtains information from mesh elements of 
similar refinement. 



Figure 8. Hole Structure Due to Level 2 Interpolation 


Level 2 interpolation is not initiated within PEGASUS until all boundary creation (hole cutting and 
outer boundary definitions) is completed. 


To implement level 2 interpolation, PEGASUS first attempts to interpolate every point in each mesh 
from all other meshes. As a result, PEGASUS performs far more point interpolations than did earlier 
versions of PEGASUS. Therefore, it is critical that the interpolation process in PEGASUS be as efficient 
and robust as possible. As in earlier versions of PEGASUS, Version 5.1 uses a Newton’s iteration 
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representation of the donor grid to provide initial guesses to the Newton's method. It was determined 
that in some cases PEGASUS was spending large amounts of computational time determining these 
initial guesses, particularly since even a sparse representation of a large grid can still represent a 
significant effort when using an exhaustive search. 


PEGASUS 5. 1 uses data structures (Alternating Digital Trees, or ADTs, Ref. 2) for each mesh to 
determine an approximate closest point in the donor mesh to the target boundary point in the recipient 
mesh. Use of the ADT is extremely fast; being a tree structure, the ADT approach can handle very large 
meshes without a significant increase in the amount of time required to find the approximate closest 
point. Once the approximate closest point is found, the Newton's method is used to iterate to the final 
stencil reference point and compute the associated interpolation coefficients. PEGASUS also uses the 
ADTs as a means to cross boundaries (e.g„ cut planes) during the Newton’s method iteration without 
requiring knowiedge of the mesh topology. 


2.5 Orphan Point Fixing 


An additional option is offered to the user that provides the ability to change second fringe "orphan" 
points back to field points. This option should only be used with flow solvers that can handle a mixture 
of single and double fringes along any interpolated boundary (e.g., OVERFLOW (Refs. 3-5) or INS3D 
(Ref. 6)). It should be noted that fixing a second fringe orphan only changes that point to a field point 
and does not affect the quality of any interpolation point that uses this second fringe orphan point as an 
interpolation stencil. Also, even if a second fringe orphan point is fixed, it is still retained as an orphan 
point within PEGASUS and efforts will be made in subsequent restarts to find a legal interpolation 
stencil for this point. 


2.6 Projection 


Problems have often arisen when two or more grids overlap at a region near a solid wall, particularly 
with the small spacing found in meshes generated to resolve boundary layers. This difficulty is due to 
the discretization of curved surfaces. Points that are on the actual surface can be found to be inside the 
discretized solid boundary of another surface. To correct this problem, mesh projection methods have 
been used successfully. Programs like PROGRD (Ref. 7) move the mesh points in one mesh relative to 
another mesh. Once the grids have been projected and stored, interpolation coefficients may be found by 
PEGASUS in the normal manner. 


PEGASUS 5. 1 has a built-in and highly automated mesh projection capability, which is based on the 
PROGRD algorithm. For a given mesh, all other overlapping meshes are projected onto its solid 
surfaces. The interpolation process then uses these projected mesh definitions to determine interpolation 
coefficients. PEGASUS does not modify the grids that will be used by the flow solver. However, each 
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definition-,. Mesii projection can be controlled through the user inputs DISTANCE, ANGLE, SHIFT, 
PINCLLDE, PEXCLLDE in the SGLOBAL and $MESH namelists. 


2.7 Lnblanking 

A feature that has been added to aid in cleaning up hole cutting is the unblanking feature. This feature is 
a manual operation that uses the SREGION/S VOLUME user inputs. Points of a specified mesh can be 
modified so as not to be hole points. This procedure is performed immediately after hole cutting. Points 
that are to be unblanked will become either a field or interpolated boundary point. 


2.8 Restarting 


In the previous version of PEGASUS, restarting was handled in a very labor-intensive manner. Each 
new restart required a new user input file that specified the changes, as well as the relationship, of the 
new added meshes to the previously added meshes. This operation was prone to errors and required a 
new approach, since the number of meshes and connections would increase prohibitively . The new r 
approach to restart requires no special knowledge of what has occurred in previous executions. The user 
is required only to modify the standard input file or the meshes, as needed. PEGASUS automatically 
determines what work needs to be performed to complete the PEGASUS execution. In this manner, 
restarting in PEGASUS is very similar to using the UNIX "make" utility. 

The procedures in PEGASUS version 5.1 have been written so create the communication data in a path 
independent manner. Path independence was not the case with restarting for PEGASUS 4. 

In order to restart a PEGASUS job from the beginning, simply remove or rename the WORK directory, 
for example: 


/bin/rm -rf WORK 
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Section 3: PROGRAM IMPLEMENTATION 


3.1 Compilation 


PEGASUS will usually be obtained in a "gzipped" tar file. When the file is unzipped and extracted, the 
top directory will contain a configure script and several subdirectories. The /peg subdirectory contains 
the source code, and the /tools subdirectory contains several PEGASUS utilities that aid in setting up the 
input, plotting, file conversion, and diagnostics. 

PEGASUS has been developed and tested on SGI and Cray platforms. Every attempt has been made to 
conform to standard Fortran 90 conventions; therefore, compilation on other platforms with compliant 
Fortran 90 compilers should be straightforward. 


To compile and install PEGASUS and all of the utility programs and scripts, follow these steps: 


1. Run the configure script. Type "configure -help" to see a list of options. Examples include: 


configure --enable-dp 
configure --enable-mips2 
configure ~-enable-mips3 
configure --enable-mips4 
configure --enable-debug 
configure --installdir=DIR 


To enable double-precision 

SGI -mips2 compiler option 

SGI -mips3 compiler option 

SGI -mips4 compiler option 

To enable debug compiler flags and output 

Executables will be installed in directory DIR 


2. Compile: type "make” 

3. Install: type "make CMD=install" 

The executable will be named pegasus5. 

It is recommended that PEGASUS be compiled and run in double precision when viscous grids are used. 
Running viscous cases in single precision has been seen to sometimes adversely affect the 
hole-generation process, and in some cases can affect interpolation as well. 


3.2 User Inputs 


L i ‘ i 

uLilily codes, which are described in .Section 3.3. which neip in the preparation or tins input. fhe details 
of these inputs are described in this section. 

The individual volume grids are stored in a sub-directory within the project directory called X_DIR. See 
Section 3.5 for a description of the directory structure used by PEGASUS5 code. The grids are to be 
stored as unformatted, single-zone, PLOT3D grid files. These grid files are named "mesh_nameA.x”, 
"mesh_nameB.x", etc, where mesh_nameA, mesh_nameB, ... correspond to the NAME variables in the 
SMESH namelists in the standard input file. 

All of the user inputs to PEGASUS are read from the standard input. All inputs are in a form similar to 
FORTRAN NAMELISTS. The following is a list of all of the NAMELISTS defined in PEGASUS . Not 
all of these NAMELISTS have to be included in the input, and the order in which the NAMELISTS 
appear in the input does not matter. Most problems will only require input to be specified in the first 
three of these NAMELISTS: 


• SGLOBAL 

• SMESH 

• SBCINP 

• SHCUT 

• SOUTER 

• SBOUNDARY 

• SSURFACE 

• $BOX 

• $ REGION 

• $ VOLUME 

• SLEVEL2 


Note that all user input record names are capitalized and all alpha-numeric names (i.e., NAME, 
ISPARTOF) are case sensitive. Examples of some these inputs is given in Appendix B. Obsolete user 
inputs (those used in version 4.x but not in version 5. 1 ) are listed in Appendix C. 

In the namelists which require grid indices, negative numbers can be used to specify indices relative to 
the last index. For example, an index of j=- 1 is equivalent to j=jmax, j=-2 is equivalent to j=jmax- 1 , etc. 

SGLOBAL 

This input sets global values for all meshes. Each record can be overridden for any individual mesh by 
using the $MESH input. However, the mesh override is limited to that individual mesh. For subsequent 
meshes, the global values are used unless specifically overridden. 


Variable 

Description 

| Default 

— 





ANGLE 


30.0 


! 


Maximum angle, in degrees, between 
the surface normal of the projected 
point and the reference surface of 
projection. If greater than this angle, 
the projection is not allowed. 


CARTX 

CARTY 

CARTZ 


Each of these inputs is a pair of 
numbers to set the min and max 
boundaries in the x, y, or z directions j 
of the Cartesian boxes used in the j 

automatic hole cutting. If they are left 
at their default values, the hole-cutting 
cartesian boundaries are set to extend 
just beyond the surfaces making up the 
hole cutter. 1 


,e30 


CENTER 


I This input can be used to limit the 
extent of the of the volume in which 
holes can be blanked out by the 
hole-cutters, and in some instances, 
can be used to help close off leaks in a 
hole-cutting surface. 


Set equal to .TRUE, to create an 
internal cell-centered grid. 


.FALSE. 


CNX, CNY, CNZ 


Integer dimensions of Cartesian grids 
(hole maps) used for hole cutting. 


I 


DISTANCE 


EPS 


EXCLBC 


A positive integer value that controls 
the maximum distance that a surface 
point can be projected. The input 
DISTANCE is the number of grid cells 
in the solid-wall normal direction used 
to compute the actual distance limit. 


Tolerance of interpolation coefficient. 
Value of interpolation coefficients less 
than 0 and greater than 1 .0 that is 
allowed for a converged stencil. 


Set equal to .TRUE, if the points 
specified as boundary conditions are to 
be excluded as interpolated boundary 
points. 


20 


0.001 


.TRUE. 
















HBFRNG 


H0LE1ST 


INCORE 


LEVEL2 


OBFRNG 


OCTLEVEL 


Set equal to 1 for single fringe or equal 


to 2 for double fringe. 


Hole boundary fringe. Set equal to 
for single fringe or equal to 2 for 
double fringe. 


If set equal to .TRUE., hole cutting, 
both manual and automatic, are 
performed prior to interpolation. 


Set equal to .TRUE, if automatic hole 
cutting is to be performed. 


Specifies if main data (e.g., meshes, 
interpolation values, etc.) are kept 
incore (=.TRUE.) or written and read 
from file (=. FALSE.) 


Set equal to .TRUE, if "level 2" 
interpolation is to be performed. 


Outer boundary fringe. Set equal to 1 
for single fringe or equal to 2 for 
double fringe. If this value is not 
specified it takes on the value of 
FRINGE. 



Integer level of octree refinements 
used to define automatic hole cutting 
boundaries. 


Set to a positive integer value, zero or 
greater. This value specifies the 
increase of the minimum hole that was 
created by an automatic hole cutter. 
For example, if OFFSET=2, all field 
points within 2 indices of the hole are 
changed to hole points. 
















ORPHF1X 

" “l 

v ,r muni so TRUE.. second fringe 
orphan points will be fixed, i.e., 
returned to field values. Certain flow 
solvers (e.g., OVERFLOW and 
INS3D) allow for mixed fringed 
conditions (single and double) and can 
use this option. 

ii 

TRUE. 

OUTER 

i 

Set equal to .TRUE, if automatic outer 
boundary specifications are to be 
performed. 

.TRUE. 

THRU 

Perform processing up to and 
including a particular operation. The 
following are the operations in the 
order that they are executed: 

projection - project meshes with solid 
boundary conditions 

octree - creates a search tree for each 
mesh 

interpolate - interpolates each 
mesh-to-mesh combination 

auto_hbound - creats automatic 
Cartesian hole cutters 

man_hboun - creates surface cutter 
based on user inputs 

auto_cut - cuts holes using automatic 
Cartesian hole cutter 

man_cut - cuts holes using manual 
hole cutting 

comp_hole - assembles composite 
holes 

spec_intl - specify points for "level 1" 
interpolation 

spec Jevel 1 - pick "best" points for 
"level 1" interpolation 

xintout 














$MESH 


Variable 

Description 

Default 

I 

QTOL 

Tolerance band for quality for 
a. This value is used when 
selecting the "best" 
interpolated boundary point 
cell-interpolation stencil pair 
(see Section 2.3 and Figure 7). 

0.0 1 or QTOL 
specified in 
SGLOBAL 

QCUTOFF 

The level of quality that is 
accepted for a. Interpolated 
boundary points with a value 
of quality below this level will 
be considered an "orphan" 
point and no interpolation 
information will be supplied. 

0.0 or the value of 
QCUTOFF in 
SGLOBAL 

HBFRNG 

Hole boundary fringe. Set 
equal to 1 for single fringe or 
equal to 2 for double fringe. 

The value of the next 
fringe value specified 
in the ascending 
order: FRINGE for 
this mesh, HBFRNG 
for SGLOBAL, 
FRINGE for 
SGLOBAL 

OBFRNG 

Outer boundary fringe. Set 
equal to 1 for single fringe or 
equal to 2 for double fringe. 

The value of the next 
fringe value specified 
in the ascending 
order: FRINGE for 
this mesh, OBFRNG 
for SGLOBAL, 
FRINGE for 
SGLOBAL 



ORPHFIX 


HANTOM 


ALFA, BETA, 
GAMA 


EPS 


FRINGE 


JINCLUDE, 

KINCLUDE, 

LINCLUDE 


Set equal to TRCH.. second 
fringe orphan points will be 
fixed, i.e.. returned to field 
values. Certain flow solvers 
le.g., OVERFLOW and 
INS3D) that allow for mixed 
fringed conditions (single and 
double) can use this option. 


Setting PH ANT OM= .TRUE . 
denotes this mesh is a 
"phantom" mesh, to be used 
only for manual hole cutting, 
and is not part of the final grid 
system. 

.FALSE. 

Euler angles (degrees) 
determining rotation of input 
meshes. ALFA, BETA, and 
GAMA are rotations about the 
Z, Y, and X coordinates, 
respectively. The coordinates 
of a mesh are transformed in 
the following sequence of 
operations: (1) scaling, (2) 
translation, and (3) rotation. 

0.0 

Tolerance of interpolation 
coefficient. Value of 
interpolation coefficients less 
than 0 and greater than 1 .0 that 
is allowed for a converged 
stencil. 

Value of EPS in 
SGLOBAL 

Set equal to 1 for single fringe 
or equal to 2 for double fringe. 

Value of FRINGE in 
SGLOBAL 

Range of J, K, and L indices, 
respectively, to be included for 
boundary points. Composed of 
two values (a starting and an 
ending value) in each of the J, 
K, and L directions, 
respectively. 

none 


















XO.YO.ZO 


OFFSET 


PROJECT 


Mesh translation factors. These 
factors are added to the original 
mesh coordinates. The 
coordinates of a mesh are 
transformed in the following 
sequence of operations: (1) 
scaling, (2) translation, and (3) 
rotation. 


Set to a positive integer value, 
zero or greater. This value 
specifies the increase of the 
minimum hole that was created 

by an automatic hole cutter. 0 

For example, if OFFSET=2, all 

field points within 2 indices of 

the hole are changed to hole 

points. 


Logical variable which 
controls whether or not 
solid-wall points in this mesh 
will actively project 
themselves onto other 



I 


solid-wall surfaces in other 


meshes. Even it this is set to 
.FALSE., this meshes solid 
walls may be used as a 
reference surfaces for other 
meshes to project onto, unless 
this mesh is listed in the 
PEXCLUDE list for all other 


Value of SGLOBAL 
PROJECT. 


meshes. 





i 

SHIFT 

1 his input is me same as the 

shift parameter in PROGRD 
and is used only if 
PROJECT=.TRUE.. If equal to 
1. move all points in J/K/L 
direction by same increment as 
surface point. If equal to 2, 
move nearest half of point in 
J/K/L direction by same 
increment as surface point. If 
equal to 3, hold opposite 
boundary fixed in J/K/L 
direction and redistribute 
interior points proportional to 
arclength distance from wall. 

I 

jy 

ANGLE 

Maximum angle between the 
surface normal of the projected 
point and the reference surface 
of projection. If greater than 
this angle, the projection is not 
allowed. 

30.0 

DISTANCE 

A positive integer value that 
controls the maximum distance 
that a surface point can be 
projected. The input 
DISTANCE is the number of 
grid cells in the solid-wall 
normal direction used to 
compute the actual distance 
limit. 

20 

PINCLUDE 

Names of meshes that are to be 
included for projection with 
this mesh. 

All meshes. 

PEXCLUDE 

Names of meshes that are to be 
excluded from projection with 
this mesh. 

No meshes. 


$BCINP 

This input is like the BCINP namelist in the OVERFLOW flow solver. This namelist specifies all of the 


beundarv condition:; for each 


PEGASUS 
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on the specification of the solid walks), and to automatically determine the outer boundary points which 
will require interpolation. It also uses these boundary conditions to determine certain mesh topologies, 
such as o-grid periodic meshes, 2D meshes, and symmetry planes. 


Variable 

Description 

Default 

ISPARTOF 

Name of MESH for which these boundary 
conditions apply. 

None 


List of boundary condition types 

I - 8 = Solid wall 

10 = O-grid periodicity (apply to first or 
last plane) 

I I = Symmetry in X 

12 = Symmetry in Y , 

13 = Symmetry in Z 

14 = Axis (J around) 

15 = Axis (K around) 

1 6 = Axis (L around) 

17 = Symmetry with no reflection plane 

2 1 = 2-D condition in Y (apply to one face) 


IBTYP 

(3 planes supplied at Y=-1,0,1) 

22 = Axisymmetric in Y, (apply to one 
face) 

(3 planes supplied, +/- 1 deg rotation about 
X) 

30 = Outflow (extrapolation) 

3 1 = Characteristic condition 

32 = Supersonic/subsonic inflow/outflow 

33 = Specified pressure outflow 

35 = Outflow (lst-order extrapolation) 

40 = Impose freestream ; 

41 = Nozzle inflow 

42 = Prescribed Q 

43 = Prescribed Q with slow start j 

43 = Impose freestream 

44 = actuator disk 

45-46 = prescribed Q/inflow-outflow 
condition 

47 = Freestream/characteristic condition 
49 = Default (no change) 

5 1 = C-grid flow-through (specify on side) 

(J is C-direction) 

52 = C-grid flow-through (specify on side) 

(K is C-direction) 

None 


IB DIR 


JBCB 


JBCE 


KBCB 


KBCE 


LBCB 


. L C-directiom 

54 = Fold-over cut flow-through (fold-over 
in J ) 

55 = Fold-over cut flow-through (fold-over 
in K) 

56 = Fold-over cut flow-through (fold-over 
in L) 

57 = C-grtd ut a wall (apply wall first) 

(J is C-direction) 

58 = C-grid at a wall (apply wall first) 

(K is C-direction) 

59 = C-grid at a wall (apply wall first) 

(L is C-direction) 

61 = Blank out region (set IBLANK=0) 

70 = Copy to 

7 1 = Copy from 

81 = Slotted wind-tunnel wall 

82 = Slotted wind-tunnel wall 
86 = Wind tunnel exit 


List of index directions: 

1 = positive J direction 

- 1 = negative J direction 

2 = positive K direction 
-2 = negative K direction 

3 = positive L direction 
-3 = negative L direction 


list of beginning j-indices 


list of beginning j-indices 


list of beginning k-indices 


list of beginning k-indices 


list of beginning 1-indices 


None 


None 


None 


None 


None 


None 


LBCE 


XSYM 


list of beginning 1-indices 


symmetry at X=0.0 


None 


1 if ibtyp=l 1 
0 otherwise 










r 
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J 

i 

symmetry at Y=0.0 

! 

! i f ibty p= ! 2 
0 otherwise 

ZSYM 

symmetry at Z-0.0 

1 if ibtvp= 1 3 
0 otherwise 

1 

X2D 

2-dimensional in X 
supply 3 planes in X=-1,0, 1 

0 

Y2D 

2-dimensional in Y 
supply 3 planes in Y=-1,0, l 

1 if ibtyp=2 1 
0 otherwise 

Z2D 

2-dimensiona! in Z 
supply 3 planes in Z=- 1,0,1 

0 


$HCUT 

This input specifies a group of meshes that be used to create a hole cutting boundary, as well as which 
meshes will be cut by this hole cutting boundary. Entries in the $BCINP namelist must be present for 
this method to work. Both solid boundary and symmetry plane information must be present in the 
$BCINP namelist for the meshes that are used for the hole cutting boundary. If no SHCUT namelist 
groups are specified and HCUT is equal to .TRUE, (in namelist group $GLOBAL), a hole cutting 
boundary will be created from all solid boundaries in all of the $BCINP namelists, and this hole cutting 
boundary will cut a hole in all meshes. There can be zero to many $HCUT namelist groups specified in 
the user inputs. 


Variable 

Description 

Default 

NAME 

Name of hole cutting group 
(must be any alphanumeric 
characters with no special 
symbols). 

none 

MEMBER 

List of meshes that make up a 
group. This group of meshes 
must define a completely 
enclosed volume. 

all meshes 

INCLUDE 

List of meshes that will have 
minimum holes cut by this 
group of meshes. 

all meshes 




II 1 

EXCLUDE 

List of meshes that \vi!! not 
have minimum holes cut bv 
this group of meshes. This 
input is only used when 
INCLUDE is not specified, i.e., 
when all meshes have been 
specified. 

1 1 

! 

none 

OCTLEVEL 

Integer level of octree 
refinements used to define 
automatic hole cutting 
boundaries. 

Global value of 
OCTLEVEL 

1 

CNX,CNY,CNZ 

Integer dimensions of 
Cartesian grids (hole maps) 
used for hole cutting. 

Global values of CNX, 
CNY, CNZ 

CARTX 

CARTY 

CARTZ 

Each of these inputs is a pair of 
numbers to set the min and 
max boundaries in the x, y, or z 
directions of the Cartesian 
boxes used in the automatic 
hole cutting. If they are left at 
their default values, the 
hole-cutting cartesian 
boundaries are set to extend 
just beyond the surfaces 
making up the hole cutter. 

Global values of 
CARTX, CARTY, 
CARTZ 


This input can be used to limit 
the extent of the of the volume 
in which holes can be blanked 
out by the hole-cutters, and in 
some instances, can be used to 
help close off leaks in a 
hole-cutting surface. 



$OUTER 

If this input is not specified and OUTER is equal to .TRUE, (in namelist group $GLOB AL) outer 
boundaries for all meshes will be specified automatically. This means that meshes which do not have 
boundary condition data given in the $BCINP namelist will have outer boundary points specified on all 
mesh faces. The $OUTER namelist group should be used only if a specific group of meshes needs to be 
included or excluded from automatic outer boundary specification. Only one $OUTER namelist group 



Variable 

Description 

Default 

" ■■ i 

INCLUDE 

List of meshes that will have outer 
boundaries determined from input 
boundary conditions. 

all meshes 

EXCLUDE 

List of meshes that will not have the outer 
boundaries determined from input 
1 boundary conditions. This input is only 
used when INCLUDE is not specified, i.e., 
when all meshes have been specified. 

none 


$BOUNDARY 


Variable 

Description 

Default 

CLOSED 

Either .TRUE, or .FALSE. Set equal to 
.FALSE, when the surfaces do not describe 
a completely closed boundary. 

.TRUE. 

ISPARTOF 

Mesh which contains boundary (must be 
any alpha-numeric characters with no 
special symbols). ISPARTOF is a string not 
exceeding ICHAR characters. 

none 

MHOLEIN 

Meshes in which the named boundary 
causes holes. MHOLEIN’s are specified as 
follows: MHOLEIN = ’mesh name 1\ 
’mesh name 2’, etc. Each mesh name in 
MHOLEIN is a string not exceeding 
ICHAR characters. 

none 

NAME 

Name of boundary (must be any 
alpha-numeric characters with no special 
symbols). A boundary name is a string not 
exceeding ICHAR characters. 

none 
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Variable 

Description 

Default 

1SPART0F 

Boundary to which surface 
belongs (must be any 
alpha-numeric characters 
with no special symbols). 
ISPARTOF is a string not 
exceeding ICHAR 
characters. 

none 

JRANGE, 

KRANGE, 

LRANGE 

Ranges of indices that 
define surface. Composed 
of two values (a starting and 
an ending value) in each of 
the J, K, and L directions, 
respectively. Maximum 
allowable ranges are: 
JRANGE: 1,JMAX 
KRANGE: 1,KMAX 
LRANGE: 1,LMAX 

Default: none 

MESH NAME 

Name of mesh that contains 
the surface (ISPARTOF) 
(must be any alpha-numeric 
characters with no special 
symbols). MESH NAME is 
used only for hole creation 
boundaries. Name is a string 
not exceeding ICHAR 
characters. 

mesh name given in 
ISPARTOF for 
SBOUNDARY. 

NVOUT 

Direction of normal 
outward from hole. NVOUT 
is required only for hole 
creation boundaries. 
Allowable values are "+J", 
"-JY+K", "-K", "+L", and 
"-L". 

none 


$BOX 





Variable 

| 

Description 

! 

Default 

ISPARTOF 

Boundary to which box belongs (must be 
any alpha-numeric characters with no 
special symbols). ISPARTOF is a string not 
exceeding ICHAR characters. 

none 

XRANGE, 

YRANGE, 

ZRANGE 

Ranges of X, Y, and Z coordinates that 
define box. Composed of two values (a 
minimum and a maximum value) in each of 
the X, Y , and Z directions, respectively. 

none 


$REGION 


Variable 

Description 

Default 

ISPARTOF 

Mesh which contains region (must be any 
alpha-numeric characters with no special 
symbols). ISPARTOF is a string not 
exceeding ICHAR characters. 

none 

NAME 

Name of region (must be any alpha-numeric 
characters with no special symbols). A 
region name is a string not exceeding 
ICHAR characters. 

none 

TYPE 

Defines type of region being specified. If 
TYPE=’HOLE\ then hole with a fringe is 
created. If TYPE=’INTR\ then an interior 
region will be blanked out but with no 
fringe. 

If TYPE=’UNBL\ specified points will not 
be hole points (This input is used to change 
hole points back to field or interpolated 
boundary points). 

HOLE 


$VOLUME 




Variable 

- ■■ — 1 

Description 

Default 

■■ 

ISPARTOF 

Region to which volume belongs (must be 
any alpha-numeric characters with no 
special symbols). ISPARTOF is a string not 
exceeding ICHAR characters. 

none 

J RANGE, 

KRANGE, 

LRANGE 

Ranges of indices that define volume. 
Composed of two values (a starting and an 
ending value) in each of the J, K, and L 
directions, respectively. Maximum 
allowable ranges are: 

JRANGE: 1,JMAX 
KRANGE: l.KMAX 
LRANGE: 1.LMAX 

none 


$LEVEL2 

If this input is not specified and LEVEL2 is equal to .TRUE, (in namelist group $GLOBAL) "level 2" 
interpolation for all meshes will be specified automatically. The $LEVEL2 namelist group should be 
used only if a specific group of meshes needs to be included or excluded from "level 2" interpolation. 
Only one $LEVEL2 namelist group (INCLUDE or EXCLUDE) can be specified. 


Variable 

Description 

Default 

INCLUDE 

List of meshes that will have "level 2 
interpolation" performed on them. 

all meshes 

EXCLUDE 

List of meshes that will not have "level 2 
interpolation" performed on them. This 
input is only used when INCLUDE is not 
specified, i.e., when all meshes have been 
specified. 

none 


3.3 Utility Codes 


Several utility codes are supplied with PEGASUS in the /tools directory: 


• peg_setup 




• peii_hoie_surt’ 

• p3d2peg 

• make_WORK 

• mesh_lister 

• clean_holes 

• peg_plot 

• peg_diag 

• peg_ho!e 

• peg_proj 

• peg_orph 

• peg_mem 

• XINtegrity 

• XIN2sp 

peg setup - This is a me nu -driven script which gives the user options to initialize a PEGASUS problem, 
remove a PEGASUS problem, prepare to restart a PEGASUS run from the beginning, or redo the 
hole-cutting from the beginning. The first option requires that the user supply a multi-zone PLOT3D 
grid file containing all of the volume grids to be used. Also, you can supply an OVERFLOW input file 
to provide the boundary conditions. It executes mesh_lister, peg_in, make_WORK, p3d2peg, and 
peg_hole_surf. This option will split up the user-supplied multi-zone grid into individual meshes and 
place them in the X_DIR directory, ready for PEGASUS5. It creates a standard input file called 
peg.in.raw to help get you started. If an OVERFLOW input file is supplied, it will add the boundary 
condition data to the input file, and will also create a PLOT3D file of the solid surfaces of the 
configuration, called hcutter.g. This file makes it easy to plot and check that the solid-wall boundaries 
form a closed surface, which is necessary for the automatic hole cutting to work. See Section 4.2 for 
more information on the use of peg_setup. 


The following six utility codes can be executed by running peg_setup, but can also be executed as 
stand-alone codes: 

peg_in [-no_add_te] - Creates a minimal PEGASUS5 input file. It is written to a file called peg.in.raw 
with the intention that it be reviewed and edited if necessary before use. It requires the meshjist file, 
which can be generated with the meshjister program. peg_in will also read the boundary condition 
data from an OVERFLOW input file and add this information into peg.in.raw The mesh names must be 
present in the OVERFLOW input file, otherwise the boundary condition information can not be used by 
PEGASUS. 

The -no_add_te flag will cause it to not add an extra solid surface (ibtyp=-l) at the trailing edge of a 
component where a flow-through boundary starts. 


peg_hole_surf - creates an unformatted, multi-zone PLOT3D grid file, containing all of the surfaces 
used in automatic hole creation. 
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make_VVORK. -creates the WORK directory and its sub-directories which are needed within the project 
directory. If these are not created before-hand, PEGASUS5 will create them when it is executed. 


meshjister - generates a list of mesh names for use by the other utility codes. The names will be 
generated from an OVERFLOW input file, if it exists; otherwise, the user will be prompted to either 
ep*er mesh names or assign default names to each grid in the user-supplied PLOT3D grid file. 


clean holes - removes the contents of /WORK sub-directories which are relevant to hole generation. 
Use this script if it is desired to re-run only the hole-generation processes again. If one one wants to 
re-run the entire PEGASUS process from the beginning, simply remove or rename the WORK directory. 


The following utility codes are supplied for diagnostic purposes: 


peg_plot [-sp] [-noorphfix] - produces an unformatted, multiple grid PLOT3D file (with IBLANKs) of 
the grid system. It reads the iblank and stencil data from the XINTOUT file, and so can only be used 
after PEGASUS5 has been run to completion (THRU = ’xintout’). Use peg_hole to plot the holes if you 
have only run with THRU = ’comp_hole’. The program prompts for the name of the output file, and for 
the level of fringes to plot, where the fringe level can be: 

1 = single fringe: blank out level-2 and higher fringe points 

2 = double fringe: blank out level-3 and higher fringe points 

3 = all: only the minimum holes are blanked out 

If the code was compiled using double-precision, one can produce a single-precision grid file by adding 
the command-line argument -sp. 

It assigns an iblank value of 101 to first level ORPHAN points so they can be plotted; for example, use 
function 3 in PLOT3D. The code will automatically assign second-level ORPHAN points an iblank 
value of 1 (a valid interior point) unless used with the -noorphfix argument, in which case the iblank will 
be assigned a value of 101. 


peg_diag - produces a diagnostic file for plotting. It can be produced after PEGASUS execution has 
been completed. The file that is created allows the user to visualize quality and cell difference parameter 
with PLOT3D or FAST (as well as any visualization software that reads PLOT3D grid and solution 
files). This file can be used with the grid file produced with peg_plot. There are no PHANTOM meshes 
included in this file. The values in this file are as follows: 


I. 

Variable 

Value or parameter 

1 

Grid point type 

Q( 1) 

Quality Parameter 
1.0 

Legal Interpolated Boundary 
Point 

Field Point 

0(2) 

Cell Difference Parameter 
0.0 

Legal Interpolated Boundary 
Point 

Field Point 

Q(3) 

Connecting Mesh Number 
Current Mesh Number 
0.0 

Legal Interpolated Boundary 
Point 

Field Point 
Orphan Point 

Q(4) 

1.0 

0.0 

-1.0 

Field Point 

Legal Interpolated Boundary 
Point 

Orphan Point 

Q(5) 

Connecting Mesh Number 
0.0 
-1.0 

Legal Interpolated Boundary 
Point 

Field Point 
Orphan Point 


pegjiole - creates a grid file with IB LANK for the hole system that has been created. PEGASUS must 
be executed through (at least) the "comp_hole" process in order to use this utility. An unformatted, 
multiple grid, with IBLANK, PLOT3D file is produced. This code must be executed in the project 
directory. 


peg_proj - creates diagnostic multi-zone PLOT3D grid and solution files, named peg_proj.g and 
peg_proj.q, which contain a zone for each solid- wall boundary surface. The grid file contains the 
grid-cooridnates at the surface points, and the q file contains fields 1, dx, dy, dz, and 1 . The dx, dy, and 
dz data are the components of the projection vector at each surface point for the projection of that point 
with the greatest magnitude. This data is best used by plotting velocity vectors at all the points, and 
coloring the vectors by velocity magnitude. Typically, it is best to scale the velocity vectors by a factor 
in the range 10 - 100. peg_proj creates a PLOT3D command file named peg_proj.com, which can be 
used to create this type of plot using PLOT3D. peg„proj also reports the maximum distance moved for 
each mesh, and onto which mesh that maximum projection occured. 


peg_orph - produces a formatted file with a list of the orphan points for all meshes or tor a specified 
mesh. The output can also be directed to the screen. 







peg_nieni - estimates memory requirements for a given PEGASUS problem. 


XINtegrity - checks the XINTOUT file produced by PEGASUS for internal consistency. This utility 
checks the accuracy of interpolations and determines whether interpolated boundary points are receiving 
information from the correct mesh elements. A PLOT3D tile is produced with [BLANKS set equal to 
the negative of the mesh number from which interpolated points receive information. Orphan points are 
identified with IBLANK=2. 


XIN2sp - converts a double precision XINTOUT file to single precision. This utility is used if it is 
desired to run the flow solver in single precision, but PEGASUS was run in double precision. 


It is recommended that all of the utility codes be executed from the project directory, as most of them 
require information from the AYORK and /X_DIR directories. 



3.4 Output File Formats 







Standard Output 


The standard output file contains the following: 


• A summary of all user inputs. 

• A final summary containing the following: 

O The total number of hole points 

O Number of interpolated boundary points for each mesh 
O Number of interpolation stencils for each mesh 
O Number of orphans in each mesh 


XINOUT 

The XINTOUT file contains all of the interpolation stencils for the fringe points and outer-boundary 
points, as well as an iblank array defining the hole-cuts. The file is unformatted, and contains the 
following four records for each mesh M. 


Record Number 

Variable Name 

1 

1BPNTS(M),IIPNTS(M),IIEPTR(M),IISPTR(M), 
MJMAX(M),MKMAX(M), MLMAX(M) 

2 

(JI(I),I=1,IINPTS(M)), 

(KI(I), I=1,IINPTS(M)), 
(LI(I), I=1,IINPTS(M)), 
(DXINT(I), I=1,IINPTS(M)), 
(DYINT(I), 1= 1 ,IINPTS(M)), 
(DZINT(I),I= 1 ,IIPNTS(M)) 

3 

(JB(I), I=1,IBNPTS(M)), 
(KB(I), I=1,IBNPTS(M)), 
(LB(I), I=1,IBNPTS(M)), 
(IBC(I),I= 1 ,IBPNTS(M)) 

4 

(IBLANK(I),I=I 1 ,12) 


where II = 1 and 12 = MJMAX(M)*MKMAX(M)*MLMAX(M) 


3.5 Directory Structure 




PEGASUS 5.1 a,,.; urn .-i’; naLUMi in a specific directory ^ructure which it will create when 

it is first executed in a particular project directory. These directories can also be created in advance using 
the make_WORK script, which is uUo executed b\ the peg_setup script. 

Execution of PEGASUS 5.1 must be performed in the project directory (see Figure 10). The user must 
supply the standard input rile, see Section 3.2 tor input descriptions) and the individual volume grids. 

The individual volume srids are stored in the X__DIR directory as unformatted, single-zone, PLOT3D 
grid files. These £nd tiles are named rnesh_narneA.x . mesh_nameB.x . etc, where mesh_nameA, 
mesh_nameB, ... correspond to the NAME variables in the SMESH namelists. 


Input: 

pegi (standard input) 
Output: 

peg jo (standard output) 

xintout 
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Figure 10. Directory Structure 
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Section 4: USING PEGASUS 


4.1 Setting Up the Input 

The following are steps that should be taken when starting a new PEGASUS problem. 

1 . Create a project directory. 

2. Prepare the volume grids as a multi-zone, unformatted, PLOT3D grid file. 

3. Optionally, prepare an OVERFLOW (Refs. 3-5) input file containing the names and the boundary 
conditions for each mesh. 

4. Run the peg_setup script. Choose option 1, which helps prepare all files necessary to run 
PEGASUS 5. You will be asked for the names of the grid and, if available, the OVERFLOW 
input file. 

5. Exit peg_setup. This will have split the volume grid file into individual volume grid files in a 
subdirectory named X_DIR, using the p3d2peg utility. This will also have generated a file named 

peg. in. raw. 

6. If an OVERFLOW input file was supplied, the peg.in.raw will contain boundary condition data in 
the BCINP namelist. If an OVERFLOW input file was not supplied, this boundary condition 
information must be added to the input file. 

7. If the boundary condition information was added by hand, run the utility peg_hole_surf to generate 
the file hcutter.g. If an OVERFLOW input file was supplied, the peg_setup script will have 
already performed this step. 

8. The hcutter.g file is a PLOT3D grid file containing the solid surfaces of the configuration. Using a 
plotting package, such as PLOT3D, carefully examine hcutter.g for gaps in the geometry that may 
affect hole generation (See Section 4.4, Operational Hints). Modify the peg.in.raw to close any 
gaps in the geometry by adding additional surfaces to the list of boundary conditions in the 
SBCINP namelists, and use IBTYP=-1 for these additional surfaces. IBTYP=-1 denotes a solid 
surface which will be used by the automatic hole cutting, but which will not be excluded from the 
list of outer boundary surface points requiring interpolation data. 

9. Once corrections have been made, rename peg.in.raw to peg.in. 

4.2 Running PEGASUS 


1. Run PEGASUS with the command: 


pegasus5 < peg.in > peg.out 




1. V. PEGASUS j.xccawv .> the XI.XTOUT i'iic will be generated, aiong with the 

PLd j.\S l S log tile. The pet;. out file will contain a table that lists, for each mesh, the number of 
interpolated points, stencils, and orphans. The log file contains information about which processes 
were executed, and the amount of execution time required. If PEGASUS was compiled in double 
precision, the XIXTOUT file will be double precision. A single precision XINTOUT tile can be 
generated by running the XIN2sp utility code. 

3. Generate a composite grid tile with the peg_plot utility. 

4. The restarting capability is a very powerful feature of PEGASUS 5. 1 . The restart function is 
modeled after the UNIX 'make" utility. A new mesh or input, or a changed mesh or input will 
cause PEGASUS to perform only the processes that are affected by these modifications. For 
example, if a new mesh is to be added to the current configuration, the mesh is simply placed in 
the /X_DIR directory and the mesh name added to the peg.in file. When PEGASUS is executed, 
the appropriate processes are performed to account for the modifications that were made. This also 
works for changes or additions in the input file peg.in. 

The restart capability allows for an incremental buildup of the PEGASUS solution. Generally, it is 
effective to initially run PEGASUS with a minimal input file, possibly simply the peg.in.raw file 
generated by peg_setup. The first run of PEGASUS will do all of the necessary interpolations. 

The interpolation step should never have to be done again, unless a mesh is changed or added. 
Subsequent runs of PEGASUS can deal with other issues, such as manipulating the hole-cutting. 

This restart capability is also very useful when PEGASUS abnormally terminates due to a system 
crash or when resource limits (e.g., CPU limits) are exceeded. Since PEGASUS retains all results 
in its work directory for all processes that have been completed, when PEGASUS is restarted, 
execution begins with the process that was executing during the termination. The user in this case 
only needs to re-execute PEGASUS; no changes in inputs are required. 


4.3 Using Externally Generated Hole Definitions 

PEGASUS uses a general purpose hole-cutting methodology that will cut good holes in the majority of 
cases for unmodified Chimera grid systems. However, there are situations where a user might want to 
use another hole-generation methodology. Or, a user might have legacy grid systems with previously 
generated hole definitions which have been successfully applied in the past. In those cases, it may not be 
necessary or desirable to use the hole-cutting methodology in PEGASUS; however, a user may still want 
to take advantage of other aspects of PEGASUS, such as level 2 interpolation. To accommodate these 
situations, a methodology has been developed whereby these externally defined hole definitions can be 
used within PEGASUS 5.1. 


To use external hole definitions, a PUOT3D file with IBLANK records must be provided, with 
IBLANK=0 specified for each hole point. The IBLANK records may be generated by any means at the 
user’s disposal. 


Run the utility code p3d2peg. The PLOT3D grid file and the mesh.list file (generated by the utility 



code meshjisten must be in the project directory. p3d2peg will u>k wa;u ;»• use the existing 

[BLANK tile for hole definition (answer in the affirmative i. li p3d2peg detects grid files in the 
/X_DIR directory, it will ask if you want to replace them. Unless one or more grids have been 
modified since the last PEGASUS run, answer "no". The 1BLANK records will be copied into the 
appropriate /WORK subdirectories. 

• Edit the peg.in file to include "HCUT=. FALSE.," in the SGLOBAL input record. 

• Run PEGASUS in the normal manner. If this is the first PEGASUS run. the entire interpolation 
process will be executed. If the interpolation process has been done earlier, the restart capability of 
PEGASUS will use the existing interpolation information to generate an XINTOUT file with the 
externally generated hole definitions. 


4.4 Operational Hints 
4.4.1 Gaps in Surfaces 


PEGASUS 5. 1 is highly automated, in the sense that much less user input is required to obtain a chimera 
interpolation file than was required in earlier versions of PEGASUS. However, automation generally 
requires making some assumptions about the input; if these assumptions are incorrect, PEGASUS may 
do something other than what the user expects or desires, resulting in varying degrees of user distress. 
One assumption is that hole cutting surfaces are watertight, in the sense that no gaps occur in the 
aggregate of the solid walls of the aerodynamic configuration. However, some common approaches to 
defining solid wall boundary conditions in flow solvers may, to PEGASUS, represent gaps which will 
adversely affect the generation of the Cartesian meshes that PEGASUS uses for cutting holes. 

Consider, for example, the solid surface generated when the OVERFLOW input used in the 3-element 
airfoil example is used to generate the peg.in file using the pegjn utility code. Figure 1 1 depicts the 
region near the corner under the back of the main wing. 



GAP IN SURFACE 



Figure 11. Surface Gaps 


In an effort to not specify the corner point twice as a solid boundary condition, the corner point is 
included in the specification of the vertical surface, but not in the horizontal surface. This is fine as far 
as the flow solver is concerned, but PEGASUS sees this as a gap in the surface. 


Figure 12 depicts the trailing edge of the wing. The upper and lower surfaces of the trailing edge are 
from two different meshes; these meshes overlap one grid point past the end of the trailing edge. Solid 
surface boundary conditions are specified up to the end of the trailing edge; beyond that point, the 
boundaries are updated via chimera interpolation. As in the corner example, PEGASUS sees a gap in the 
surface. 


The effects of these gaps become apparent during the hole-cutting process. Recall from Section 2. 1 that 
the outside region is defined by marching inward from the corners of the Cartesian mesh; any 
unidentified Cartesian element that is adjacent to an outside element must also be an outside element. It 
is clear that if the Cartesian mesh is fine enough, outside elements can "leak" through the gaps in the 
surface and misidentify elements inside the body as outside elements. Therefore, no hole will be cut. 



OPENING AT TRAILING EDGE 



Figure 12. Trailing Edge Opening 


The remedy for gaps in the surfaces depends on the nature of the openings. In the case of the corner gap, 
the peg.in file can simply be modified to extend the solid boundary of the horizontal surface to the 
corner. The trailing edge gap, however, requires the addition of another boundary surface in the $BCINP 
namelist, and use IBTYP=-1 for this additional surface. IBTYP=-1 denotes a solid surface which will be 
used by the automatic hole cutting, but which will not be excluded from the list of outer boundary 
surface points requiring interpolation data. 

There is always a possibility that a small gap, such as that shown in the airfoil trailing edge in Figure 12. 
If the gap is smaller than the vertical spacing of the Cartesian mesh, the surface is closed as far as the 
Cartesian mesh is concerned. In some cases, the gap could be up to twice the corresponding spacing of 
the Cartesian mesh and still result in a good hole definition. 

Other methods for fixing surfaces with gaps is with the use of auxiliary (i.e., "phantom") meshes. These 
meshes are designated in the standard input file in the $MESH namelist by setting PHANTOM=.TRUE. 
The solid surfaces of a PHANTOM mesh are to be designated in the BCINP namelist for the phantom 
mesh. The phantom meshes do not participate in the interpolation, and do not appear in the output 
XINTOUT file. 

Finally, problems with automatic minimum hole cutting can be remedied with the use of the manual 
hole cuts, which use the hole cutting algorithm in PEGSUS 4.1. Manual hole cutting surfaces are 
specified using the $BOUNDARY and $SURFACE namelists. 


4.4.2 Thin Trailing Edges 



The hole cutting algorithm is a very general approach, and can accommodate unmodified chimera grids 
systems with overlapped collar grids and multiple interior regions (which occur, for example, in the 
3-element airfoil example). However, modern chimera systems often have a very large size range of 
surface features that the Cartesian hole cutting meshes must resolve. In the case of very fine features on 
a large configuration, the resolution required may be excessive. Consider, for example, the trailing edge 
of slat in the 3-element-airfoil example. Figure 13 depicts a point inside the 



Figure 13. Hole Point Identification with Thin Trailing Edge 


trailing edge of the slat. The solid surfaces of the trailing edge are encompassed by fringe elements, only 
one of which is shown for purposes of clarity. Recall from Section 2.1 that points inside fringe elements 
use a "line-of-sight" algorithm to determine their status. To reiterate, if a clear line-of-sight exists from a 
point in a fringe element to another non-fringe element, then the point adopts the status of the non-fringe 
element (either "outside" or "inside"). However, the point in this example may not have a clear 
"line-of-sight" to a non-fringe element, and therefore cannot be identified in this manner. If a point 
cannot be identified, it is assumed to be an outside (i.e., field) point. If the field point is next to a hole 
point, it will become a candidate for interpolation; being inside the slat, there will be no viable 
interpolation stencil, and the point will be an orphan. 

Fortunately, points in thin trailing edges can often be blanked as desired using the OFFSET input in the 
SGLOBAL or $MESH input records. Setting OFFSET = 2, for example, grows hole regions by blanking 
any non-hole point within 2 indices of a hole point. In this manner points in thin trailing edges that are 
destined to be orphaned may be correctly blanked. 

Also, multiple hole cutting groups can be used in PEGASUS. Separate hole cutters can be designated, 
for example, for the slat, wing, and flap separately. The Cartesian cutters are automatically sized 




automatically increase me Jeiauil resolution oi uie cartesian noie-cutting mesh in the region of the 
trailing edge. An example of input which invokes multiple hole cutters is reproduced in Appendix B. 
Finally, the default logical dimensions of the Cartesian hole cutter may be increased until there is 
sufficient mesh resolution to resolve the trailing edge. 


4.4.3 Projection 

The value for DISTANCE in the SGLOBAL and SMESH namelists will control how many points can 
be projected. If DISTANCE is too large, it is possible for PEGASUS to project points too far, and 
actually move them from one component to another. It is advisable to routinely check the results of the 
projection by checking the maximum distance projected. This is reported in the log file, and is also easy 
to check by running the utility peg_proj. 
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octane13:rogers 

stdin 

Mon Aug 28 07:01:33 2000 
hpl 1 / hp8apple 




c u i_ uii e — -3 : rogers 

IjOO i 

stdin 

hpll 

octanel3 : rogers 

Job: 

stdin 

hpll 

octanel3 : rogers 

Job: 

stdin 


Date: Mon Aug 28 07:01:33 2000 
Date: Mon Aug 28 07:01:33 2000 
Date: Mon Aug 28 07:01:33 2000 


hpll octane!3 : rogers 


Job: stdin Date: Mon Aug 28 07:01:33 2000 
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Section 6: APPENDICES 


APPENDIX A: Glossary 



Donor Mesh 


Field Point 


Fringe Point 


Hole Boundary 


Hole Boundary Point 


Hole Map 


Hole Point 



Definition 


Mesh that provides information by interpolation 
to another mesh. 


A point within computational 
domain(IBLANK=l) that is updated by the flow 
solver or boundary conditions. 


A hole boundary point. 


Collection of hole points, i.e., points in a mesh 
that have been blanked. 


Collection of hole boundary points (fringe points) 
surrounding a hole. 


Point on a hole boundary, by definition a point 
between a hole point and field point. 


A cartesian grid that has been developed to cut 
minimum holes. 


A mesh point that has been specified to be within 
a hole. This point is considered to be out of the 
computational domain (IBLANK=0). 

















Example o l PKGa.>LS input Pile ipeg.in) 
(3-Element Airfoil) 


EXAMPLE 1 - Automatic Hole Cutting and Outer Boundary Specification: 

The following input file relies on the automatic definition of the minimum hole cutter, which uses all of 
the solid walls to cut holes out of all of the zones. 


SGLOBAL 
PROJECT = ,F„ 

LEVEL2 = T.. 

FRINGE = 2. 

CNY = 3, 

OFFSET = 2, 

SEND 
SMESH 
NAME = ’slat’, 

SEND 

SMESH 

NAME = ’main’, 

SEND 

SMESH 

NAME = 'flap 1\ 

SEND 

SMESH 

NAME = ’flaptel’, 

SEND 

SMESH 

NAME = ’wake’, 

SEND 

SBCINP 

ISPARTOF= ’slat’, 

IBTYP = 5, 51,21, 

IB DIR = 2, 2, 3, 

JBCS= 15, 2, 1, 

JBCE = 107, 14,-1, 

KBCS = 1, 1, 1, 

KBCE = 1, 1,-1, 

LBCS = 1,2, 1, 

LBCE = -1, -2, 1, 

SEND 

SBCINP 

ISPARTOF= ’main’, 

IBTYP = 5, 47, 47,21,47,-1, 
IBDIR = 2, 1,-1, 3,-2, 2, 

JBCS = 121, 1, 321, 1, 1, 231, 
JBCE = 231, 1, 321, -1, -l, 232, 
KBCS = 1, 1, 1, 1,-1, 1, 

KBCE = 1, 101, 101, -1,-1, 1, 
LBCS = 1,2.2, 1,2, 1, 

LBCE = -1,-2, -2, 1,-2, -1, 



>FND 

sD.c 

ISPAR TOP = 'Hap!'. 
[BTYP = 5, 21, 

[BDIR = 2, 3. 

JBCS = 21. I, 

JBCE = 121. -1. 

KBCS =1.1. 

KBCE = l.-l. 

LBCS =1.1. 

LBCE = - 1 . I. 

SEND 

SBCINP 

[SPARTOF = ’flaplte', 
[BTYP = 5, 21. 

IBDiR = 1.3, 

JBCS = 1, 1, 

JBCE = 1,-1, 

KBCS = 4. 1. 

KBCE = 28. - !. 

LBCS = 1, 1, 

LBCE = - 1 , 1, 

SEND 

SBCINP 

ISPARTOF = ’wake’, 
IBTYP = 5, 5,47, 21, -L 
IBDIR = 1,-2. -1,3, -2, 
JBCS = 1, 1, 151, 1, 61, 
JBCE - 1,61, 151,-1,62, 
KBCS = 1.61, 1. l.-l. 
KBCE= 61, 61, 61, -1, -1, 
LBCS = 1, 1,2, 1, 1, 
LBCE = -1,-1, -2, 1,-1, 
SEND 


EXAMPLE 2 - Automatic Hole Cutting Using Three Separate Hole Cutters: 

By adding these $HCUT namelist inputs to the input file in the previous example. These namelists 
define three hole-cutting surfaces, one for the slat, one for the main wing, and one for the flap. 


SHCUT NAME = ’slathole’, 
MEMBER = ’slat’, 

INCLUDE = ’main’, 

SEND 

SHCUT NAME = ’ winghole’, 
MEMBER = main’, ’wake’, 
INCLUDE = ’slat’, ’ flap 1 ’ , 
SEND 

SHCUT NAME = ’flaplhole’, 
MEMBER = ’flap 1’, ’ flaplte’, 
INCLUDE = ’main’, ’wake’, 
SEND 



Obsolete Inputs 


The following inputs were available in versions of PEGSUS 4. 1 but are no longer accepted as inputs 
PEGASUS 5.1. 


NAMELIST 

Obsolete Input 

$GLOBAL 

QUALITY 

$MESH 

LINK, ADDLINK. ADDHOLE, QUALITY 
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